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ABSTRACT 



We present an application of the da Cunha, Chariot & Elbaz (2008) model of the spectral energy distribution (SEDs) of galaxies from 
the ultraviolet to far-infrared, to a small pilot sample of purely star-forming Ultra-Luminous Infrared Galaxies (ULIRGs). We interpret 
the observed SEDs of 16 ULIRGs using this physically-motivated model which accounts for the emission of stellar populations from 
the ultraviolet to the near-infrared and for the attenuation by dust in two components: an optically-thick starburst component and the 
diffuse ISM. The infrared emission is computed by assuming that all the energy absorbed by dust in these components is re-radiated at 
mid- and far-infrared wavelengths. This model allows us to derive statistically physical properties including star formation rates, stellar 
masses, as well as temperatures and masses of different dust components and plausible star formation histories. We find that, although 
the ultraviolet to near-infrared emission represents only a small fraction of the total power radiated by ULIRGs, observations in this 
wavelength range are important to understand the properties of the stellar populations and dust attenuation in the diffuse ISM of these 
galaxies. Furthermore, our analysis indicates that the use of mid-infrared spectroscopy from the Infrared Spectrograph on the Spitzer 
Space Telescope is crucial to obtain realistic estimates of the extinction to the central energy source, mainly via the depth of the 9.7 '-pm 
silicate feature, and thus accurately constrain the total energy balance. Our findings are consistent with the notion that, in the local 
Universe, the physical properties of ULIRGs are fundamentally different from those of galaxies with lower infrared luminosities and 
that local ULIRGs are the result of merger-induced starbursts. While these are well-established ideas, we demonstrate the usefulness 
of our SED modelling in deriving relevant physical parameters which provide clues to the star formation mode of galaxies. 

Key words. Galaxies: evolution - galaxies: fundamental parameters - galaxies: starburst - galaxies: ISM. 



1. Introduction 

Ultra-Luminous Infrared Galaxies (ULIRGs) are galaxies with 
infrared luminosities higher than 10 12 Lq. They emit most of 
their energy (over 90 per cent) in the infrared, which implies 
that they are heavily dust-obscured. To power their huge infrared 
luminosities, these galaxies must be undergoing intense star for- 
mation, forming new stars at rates of the order of 100 Mq yr _1 . 
In the local Universe, such large starbursts are thought to re- 
sult from major interactions or mergers of gas-ric h galaxies 
(e.g. iMurphv et all 120011: ISanders & Mirabell Il996l and refer- 
ences therein). This is supported by the fact tha t local ULIRGs 
have typically disturbed morphologies (e.g . IClements et al.l 
Il996t iDuc et al.lfl997b iRigopoulou et al] 1 19991) . Another possi- 
ble source of energy in ULIRGs is the accretion of large quan- 
tities of gas onto a central super massive black hole in an ac- 
tive galactic nucleus (AGN; e.g. ISanders et al.lll988l) . In such 
cases, dust in an optically-thick torus surrounding the AGN 
absorbs large amounts of energy and re-radiates it at longer 
wavelengths, contributing to the large infrared luminosities of 
ULIRGs. In most cases, there is likely to be a combination of 
both star formation and AGN contributing to the emission from 
ULIRGs, and the relative contribution by the two processes has 



been extensively investigated (e.g. Lutz et all l998t Genze et all 
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1998; Laurent et al. , 2000; Armus^t_alJ|2007; Spoon et al. 2007; 
IVeilleux et al.ll2009l) . 

Understanding the formation and evolution of ULIRGs is 
key to understanding the cosmic evolution of galaxies. ULIRGs 
play a significant role in galaxy evolution at high redshifts and 
they are major connibutors to the infrared lumi nosity and star 
formation density at redshift s z ~ 1 — 2 (e.g. ICharv & Elbazl 
120011: iLeFloc'h etall 120051) . Local ULIRGs may also pro- 
vide analogs to dusty galaxy populations at high red shifts 
(e.g. iBlain etaDl2002t Daddi et alJl2007t iDev et al.ll2008l) . To 
understand the evolution of ULIRGs in a cosmological context, 
the study of the physical properties of these system both at low 
and high redshifts is required. 

With the advent of a number of sensitive space observa- 
tories and new instrumentation on ground-based telescopes, a 
wealth of new data on local samples of ULIRGs have been com- 
piled over recent years. The analysis of such data has focused 
mainly on the infrared range, where ULIRGs emit most of their 
energy. Several models have been proposed to inter pret the in- 
frared spectral energy distributions of ULIRGs (e.g. [ Klaas et all 
12001b iFarrah et alJ2003t IVega et all2008t Marshall e t al.ll2007l) . 
However, such models are not easily applicable in the statisti- 
cal study of large samples of galaxies and they often neglect the 
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optical emission of ULIRGs, which may carry important infor- 
mation about their intermediate-age and old stellar populations 
and hence the ir star formation histori es. 

Recently, Ida Cunha et al.l d2008l) have presented a simple, 
physically motivated model that allows the interpretation of the 
integrated ultraviolet, optical and infrared spectral energy dis- 
tributions (SEDs) of galaxies in te rms of fundamental ph ysical 
parameters. This model was used in lda Cunha et alj j2008) to de- 
rive median-likelihood estimates of parameters such as the star 
formation rate, stellar mass, dust attenuation and dust mass from 
the observed SEDs of 66 galaxies in the Spitze r Infrared Nearby 
Galaxy Survey ( SING S, iKennicutt et alj2003l) . In another study, 
Ida Cunha etatl(l2010h have used this model to investigate the 
relation between star formation and dust content of a sample of 
3258 low-redshift IRAS galaxies detected by Sloan Digital Sky 
Survey (SDSS) with complementary observations by GALEX 
(ultraviolet) and 2MASS (near-infrared). These studies have pro- 
vided valuable insight into the physical properties of local galax- 
ies with moderate tot al infrared luminositi es (Lir < 10 12 L e ). 

The method of Ida Cunha et al. was designed to in- 

terpret SEDs with sufficiently wide wavelength coverage (from 
the UV to the far-IR) to allow for accurate extinction corrections 
and energy balance. The main assumptions are that the energy 
input originates only from stars and all energy absorbed by dust 
in the ultraviolet and optical is re-radiate d in the infrared . In this 
paper, we investigate how the model of Ida Cunha et al.l (2008) 
can be calibrated to interpret multi-wavelength observations of 
star-forming ULIRGs, which have much higher infrared lumi- 
nosities and for which only a small fraction of the energy is de- 
tected in the ultraviolet and optical. This is the first step of a 
larger study in which we aim at expanding our methodology to 
interpret the physical properties of samples of galaxies covering 
a wide range of infrared luminosities, star formation and AGN 
activities, such as the Great O bservatories All-Sky LIRG Survey 
(GOALS: lArmus et all 120091) and the local ULIRG sample of 
iDesai et al.l d2007l) . 

In the present study, we start by analysing the SEDs of 
16 purely star-for ming ULIRGs from the primary sample of 
IDesai et al.l (120071) . In a forthcoming paper, we will extend our 
analysis to galaxies with AGN contribution. 

This paper is organized as follows. In Section|2] we describe 
our sample of local ULIRGs and the multi-wavelength data used 
in our study. We detail our method to derive statistical constraints 
on the physical parameters of our galaxies from fits to their ob- 
served spectral energy distributions in Section [3] In Section |4] 
we provide the fits to the SEDs of the galaxies and the estimates 
of their physical parameters. In Section [5] we show that our re- 
sults are consistent with previous findings regarding the forma- 
tion and evolution of ULIRGs, and compare the physical prop- 
erties of local galaxies in two bins of total infrared luminosity: 
10 < log(L IR /L Q ) < 12 and 12 < log(L IR /L Q ) < 12.5. Our sum- 
mary and conclusions are presented in Section [6] When neces- 
sary, we adopt the following cosmology: Hq = 70 km s _1 Mpc -1 , 
Q M = 0.30 and Q A = 0.70. 

2. The sample 

We select our galaxies from a large sample of 107 local 
ULIRGs with 0.02 < z < 0. 93 observed in the 5 - 38 ami 
regio n with the Spitzer/IRS (lArmus et alj 120071: IDesai et al] 
2007). The s ources were primarily selected from the IRAS 
1 Jy survey dKim & Sanders! Il998l) . the IRAS 2 Jy survey 
dStrauss et alj 19921) and t he FIRST///M5' radio-far-infrared sam- 
ple of Stanford et al] d2000h . We note that preliminary theoretical 
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Fig. 1. Redshift distribution of the 16 star-forming ULIRGs used 
in this study (red histogram). For comparison, the grey histogram 
shows the r edshift distribution of the primary sample of 107 local 
ULIRGs of lDesai et al.l d2007l) . 



interpretat ion of infrared obs ervations of this sample was pre- 
(120071) using the spectral decomposition 



formed by Armus et al 
method of Mars hall et al.1 (120071) 



2.1. Selection of our star-forming ULIRGs 

Given the theoretical tools we have in hand, in the present study, 
our goal is to explore how the physical properties of purely 
sta r-forming ULIRGs ca n be characterized using the model 
of Ida Cunha et all d2008l) . Modelling the SEDs of ULIRGs 
presents a challenge since theses systems have infrared lumi- 
nosities a factor of 10 greate r than what had been studied so far 
dda Cunha et al.l I2008L l2010h and very high infrared-to-optical 
ratios. Typically more than 90% of the energy output of ULIRGs 
is emitted in the infrared. This energy originates not only from 
young massive, heavily obscured starbursts, but in some cases 
from an AGN. Identifying the presence of an AGN in the dust- 
enshrouded nucleus of a ULIRG is often challenging, in partic- 
ular in the optical, where the extinction effects are more preva- 
lent. The best diagnostic is the detection of hard X-ray emis- 
sion, or high ionization lines which cannot be excited by stel- 
lar photospheric emission. The task is somewhat easier in the 
infrared, and various methods have been developed to reveal 
an AGN and quantify i ts contribution to the bolometric emis- 
sion of galaxy (see e . g.|Genzel et al.lfl 998: Laurent et al. 2000; 
l lmanishi et aLj l2006t lArmus et al.l 12001 Icharmandarisl 12008: 
IVeilleux et alj |2009l). In the present study, we focus on the mod- 
elling of a sample of ULIRGs with no indication of significant 
AGN contribution. The inclusion of galaxies with AGN activity 
will be the subject of a forthcoming study. 

We select the star-f orming ULIRGs to be studied here from 
the primary sample of Desa i et al.l d2007l) . Our selection pro- 
cess of the 16 systems is not designed to be complete. Following 
the motivations described above we have identified galaxies 
for which none of the well-known optical, near-infrared and 
mid-infrared spectroscopic diagnostics show indication of an 
AGN dominating in the optical or mid-infrared. Our 16 galax- 
ies have strong polycyclic aromatic hydrocarbon (PAH) emis- 
sion, in particular at 6.2 /zm, and no [NeV]/U4.33^m emission, 
or a mid-IR continuum indicative of an AGN. Additionally, all 
the ULIRGs in our sample can be classified as 'cool' according 
to their IRAS 25- to 60-^m flux ratios (i.e. F25/F60 < 0.2, see 
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Table [TJ, which has been taken as an indication of no significant 
AGN contribution to the mid-IR emission (e.g. ISanders et alJ 
Il988t iFarrah et al.l 120071) . We note that even though two sys- 
tems (IRAS 121 12+0305 and IRAS 16334+4630) are classified 
as LINERs from optical s pectroscopy in the literature dKim et ail 
1998i lVeilleu~ et al.ll2002l) . this does not necessarily imply the 
presence o f an AGN but rather evidence o f starburst-driven 
shocks (see lKim et alJIl 9981 : iLutz et al.lll999l) . Furthermore, de- 
tailed analysis of the infrared emi ssion of these galaxies using 
the method o f iMarshall etaP d2007l) shows that the contribution 
of an AGN to their infrared luminosity would be negligible (in 
the case of IRA S121 12+0305, this is further supported by the 
SED analysis of IFarrah et al.ll2003l who estimate that an AGN 
would contribute by at most 0.1% to the total infrared luminos- 
ity of this galaxy). 

In Fig. [T| we compare t he re dshift distribution of the pri- 
mary sample of iDesai et "all (l2007h with that of our sub-sample 
of 16 star-forming ULIRGs. We note that our galaxies are at typ- 
ically lower redshift than the bulk of the galaxies from the pri- 
mary sample due to the availability of higher signal-to-noise data 
for nearby systems which allowed us to select them as purely 
star-forming. The basic properties of our ULIRGs are listed in 
Table □ 




A/jUm (rest-frame) 

Fig. 2 . IRS spectra of the 16 ULIRGs analyz ed in this work (grey 
lines: iDesai et alJfeOOTt lArmus et al. 2007). These mid-infrared 
spectra show strong PAH features in the range 6 to 17 fan, as 
well as a rising continuum from hot dust out to 35 pm, indicat- 
ing strong star formation in these systems. The strong absorp- 
tion feature at A ^ 10 pm is due to silicates and indicates high 
dust optical depths in these galaxies. For comparison, we plot the 
starburst template of Brandl et al. (2006) in red. All the spectra 
are normalized to unity at 15 pm. The strengths of the silicate 
absorption feature and the 6.2-yimi PAH equivalent widths of all 
spectra are presented in Table [2] 



2.2. Mid-infrared Spitzer/IRS spectroscopy 

All ULIRGs in the sample of Des ai et al ] (12007b were observed 
in staring mode with both subslits of the Short-Low (SL) and 
Long-Low (LL) modules of the IRS. The final spectra have a 
spectral resolution typically between 60 and 120 over the wave- 
length range from 5 to 38 /mi. More de tail on the IRS obs erva- 
tions and data reduction can be found in Des ai et al. I J20071) . 

In Fig. |2] we show the mid-infrared IRS spectra of the 16 
ULIRGs studied in this paper. The well-known PAH emission 
features, as well as silicate absorption, molecular hydrogen lines 
and neon fine structure lines are clearly visible. It is interesting to 
compare the mid-infrared spectra of our ULIRGs with those of 
typical local starbursts of lower infrared luminosit y. Therefore, 
we al so plot the average starburst IRS spectrum of B randl et"a"D 
(2006) in Fig. [2] In general, the shapes of the mid-IR emission 
from our star-forming ULIRGs and the starburst template are 
remarkably similar, except for two striking differences: (i) most 
of our ULIRGs present much higher silicate absorption, meaning 
higher dust optical depths; (ii) the rising continuum from hot 
dust at A <: 20 pm is steeper for the U LIRGs, indicating that 
dust is heated to higher temperatures (see lDesai et al ] |2007T for a 
more detailed analysis of the mid-infrared spectra of ULIRGs). 

Given the small spatial scales in which star for mation takes 
place , our ULIRGs are essentially unresolved (e.g. lSoifer et al.l 
I2001L Diaz-Santos et al. 2010), and therefore our IRS spectra 
are representative of the total emission of the galaxies, and ca be 
used to characterize the global physical properties of our galax- 
ies. We note that we do not aim at fitting the detailed shape of the 
mid-infrared spectra with the model used in this paper, but rather 
the global SEDs from the ultraviolet to the far-infrared. We use 
the IRS spectra to help select star-forming galaxies using mid-IR 
AGN diagnostics, and also to extract essential information such 
as the depth of the silicate feature at 9.7 //m (Section [3Tj2| and 
mid-infrared broad-band colours (Section 12.3b to provide gen- 
eral information about the shape of the mid-infrared continuum 
of the galaxies. 



2.3. Ancillary multi-wavelength photometry 

We compile available photometric data in the ultraviolet, op- 
tical and infrared in order to obtain a sampling of the ob- 
served spectral energy distribution of each galaxy as com- 
plete as p ossible, which i s cruci al to fully exploit the poten- 
tial of the da Cunha et al. (2008) model. We include observa- 
tions obtained with GALEX, SDSS, 2MASS, HST, ISO, IRAS 
and Spitzer found in the NASA/IPAC Extragalactic Database 
(NED). Optical fluxes for IRAS 17208-0014, IRAS 19297-0406 
and IRAS22491-1808 were kindly provided ahead of publica- 
tion by D. Sanders and V. U. For consistency, all the selected 
fluxes correspond to the total flux from the galaxies in each band. 
Additionally, we compute mid-infrared fluxes in the IRAS \2-pm 
and 25-fim, ISO 15-pm, Spitzer/IRAC 5.8-/im and 8.0-pm, and 
MIPS fluxes 24-yum bands by convolving the observed IRS spec- 
tra with the respective filter response functions. 

Several of our ULIRGs appear do uble in the optical and 
near-IR imaging (e.g. lVeilleux et al.l2002l) . However, in all cases 
we have Spitzer IRAC and/or IRS \6-pm images which allow 
us to identify multiple mid-infrared emission components more 
clearly. These images reveal that, in the case of double sys- 
tems, either the two components are too close and both are in- 
cluded within the 3.6 arcsec width of the IRS slits, or one of the 
two clearly dominates, contributing with over 90% of the mid- 
infrared flux. We have checked that in these cases the IRS spectra 
are taken on the mid-infrared dominant source, which continues 
to remain unresolved at longer infrared wavelengths. In addi- 
tion, for consistency, we use ultraviolet to near-infrared fluxes 
extracted using large apertures which include the total flux from 
both sources in the case of double systems. 

In Appendix [A] we list all the fluxes used in this work from 
the ultraviolet to the far-infrared. 
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Table 1. Basic properties of the 16 ULIRGs analyzed in this paper, (a) IRAS name; (b) right ascension (from NED); (c) declination 
(from NED); (d) redshift; (e) total infrar ed luminosity (8 - 1000 Ai m), computed using the IRAS fluxes at 12, 25, 60 and 100 /im 
(Table lAJl l, following the prescription in lSanders & Mirabel! (I 1996b : (f) ratio of the IRAS fluxes at 25 and 60 /im. 



Galaxy 


RA. (J2000) 


Decl. (J2000) 


z 


log(L IR /L ) 


Fis 1 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


IRAS00 199-7426 


00 h 22 m 07 s .0 


-74°09'42" 


0.0963 


12.29 


0.078 


IRAS0 1494- 1845 


01 h 51 m 51 s .4 


-18°30'46" 


0.1579 


12.20 


0.059 


IRASZ02376-0054 


02 h 40 m 08 s .6 


-00°42'04" 


0.4097 


12.39 


0.024 


IRAS041 14-51 17 


04 h 12 m 44 s .2 


-51°09'41" 


0.1245 


12.18 


0.035 


IRAS06009-7716 


05 h 58 m 37 s .O 


-77°16'39" 


0.1169 


12.00 


0.083 


IRASF08208+321 1 


08 23 54 .6 


+32 02 12 


0.3957 


12.32 


0.024 


IRASF10156+3705 


10 h 18 m 34 s .5 


+36°49'52" 


0.4897 


12.46 


0.016 


IRAS 10565+2448 


10 h 59 m 18 s .l 


+24°32'35" 


0.0431 


12.03 


0.105 


IRAS 121 12+0305 


12 h 13 m 46 s .O 


+02°48'38" 


0.0726 


12.27 


0.057 


IRAS13120-5453 


13 h 15 m 06 s .3 


-55°09'23" 


0.0307 


12.24 


0.073 


IRAS 16334+4630 


16 h 34 m 52 s .6 


+46°24'53" 


0.1909 


12.37 


0.067 


IRAS 17208-0014 


17 h 23 m 21 s .9 


-00°17'01" 


0.0430 


12.39 


0.050 


IRAS 19297-0406 


19 h 32 m 2P.2 


-03°59'56" 


0.0857 


12.38 


0.071 


IRAS 19458+0944 


19 h 48 m 15 s .7 


+09°52'05" 


0.0999 


12.31 


0.064 


IRAS20414-1651 


20 h 44 m 18 s .2 


-16°40'16" 


0.0870 


12.16 


0.061 


IRAS22491-1808 


22 h 51 m 49 s .2 


-17°52'23" 


0.0772 


12.15 


0.097 



3. Method 

3.1. Spectral energy distribution modelling 

We use the simp le, physically motivated model of 
Ida Cunha et alJ (I2008I) to interpret the mid- and far-infrared 
spectral energy distributions of galaxies consistently with the 
emission at ultraviolet, optical and near-infrared wavelengths. 
This model was previously calibrated using samples of star- 
forming galaxies with moderate star formation rates and infrared 
luminosities. Here we investigate how well this prescription 
can account for the SEDs of local star-forming ULIRGs, which 
have much higher infrared-to-optical ratios. In this section, we 
briefly recall the main features of the model of Ida Cunha et alJ 
(2008), and describe the modifications we have performed to 
the original model in order to make it more suitable for the 
interpretation of the multi-wavelength emission by ULIRGs. 

3.1 .1 . Stellar emission 

We compute the emission by stars i n gal axies using the lat- 
est version of the lBruzual & Charlod (I2003) population synthe- 
sis code (Chariot & Bruzual, in preparation). This code pre- 
dicts the spectral evolution of stellar populations in galaxies 
from far-ultraviolet to far-infrared wavelengths and at ages be- 
tween 1 x 10 5 and 2 x 10 10 yr, for different metallicities, initial 
mass functions (IMFs) and star formation histories. We adopt the 
IChabrierl (120031) Galactic-disc IMF. 

3.1.2. Dust attenuation 

The ultraviolet, optical and near-infrared emission from stars is 
attenuated by dust. We ac c ount f or this using the two-component 
model of lCharlot&Falll (l2000l) . This model includes the fact 
that stars are born in dense molecular clouds with finite life- 
times fflc; after that time, these birth clouds dissipate as a result 
of the pressure caused mainly by strong stellar winds and su- 
pernovae explosions, and stars migrate to the ambient (diffuse) 
ISM. Thus, the light produced by stars younger than tsc is atten- 



uated by dust in the birth clouds and in the ambient ISM, while 
the light produ ced by older stars is attenuated only by dust in the 
ambient ISM. ICharlot & Fahl (120001) have calibrated this model 
using a sample of local, ultraviolet-selected starburst galaxies of 
low infrared luminosities, and show that a birth cloud lifetime of 
fee = 10 7 yr can account for the observed line and continuum 
emission properties of normal star-forming galaxies. 

Observations of local ULIRGs show that these galaxies are 
forming stars at high rates in a very concentrated central region 
with typical scal e of at most a few kpc (e . g. CO interferometer 
obser vations by iBrvant & Scovillj 1 19961 iDownes & . Solomon! 
1998; also observations of Arp220 bv lSakamoto et al.ll2008l) . In 
this scenario, the central region of the ULIRG would be sim- 
ilar to a huge, optically-thick molecular cloud wi th a lifetime 
typica lly larger than 10 7 yrs. Further evidence by Lahuis et al. 
(2007) suggests the presence of high abundances of warm, 
dense gas, associated with deeply embedded star formation in 
ULIRGs. In this environment, Hn regions are prevented from 
expanding by large pressure gradients of gravity, thus increas- 
ing the lifetime of the star formation pr ocess. Therefore, in th e 
framework of the attenuation model of ICharlot & Fa ll (2000), 
we adopt tsc = 10 8 yr for ULIRGs. The interstellar medium 
is still described using two components. The first is the 'birth 
cloud' component, which accounts for the central, optically- 
thick starburst powering most of the infrared emission from 
these galaxies. The second component, the 'diffuse ISM', has 
smaller effective optical depths and is heated by stars older than 
tsc- This component is necessary to account for the small but still 
present ultraviolet and optical emission from ULIRGs, which 
cannot arise from the central starburst due to the extremely high 
optical de pths (see also spectro scopic evidence supporting this 
scenario in lSoto & MartinboioT) . 

We use an 'effective absorption' curve for each component, 
f A oc A~ n , with the slope n reflecting both the op tical properties 
and th e spatial distribution of the dust. Following lCharlot & Falll 
(l2000h . we adopt for the diffuse ISM 

f] SM = ^f v (^/5500 Ar - 7 , (1) 
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where fy is the total effective V-band absorption optical depth 
of the dust seen by young stars inside birth clouds, and fi = 



iISM 



/(? 



BC 



+ T 



ISM 



) is the fraction of this contributed by dust in the 



surrounding diffuse ISM. For the stellar birth clouds, we adopt: 



ff c = (1 -//)fv01/5500A) 



-1.3 



(2) 



The slopes of the effective absorption curve in the birth 
clouds and in the diffuse ISM are different in order to account 
for the different spatial distributions of dust. For the birth cloud 
component, we adopt n = 1.3, which corresponds to a fore- 
ground dust screen. This should be adequate to characterize 
the attenuation of optically-thick dust in the concentrated dense 
central starburst region of ULIRGs. The diffuse ISM is bet- 
ter described by a ran dom distribution of discrete dust clouds 
dCharlot & Falll I2000I) . which can be accounted for adopting 
n = 0.7 (i.e. a 'greyer' effective absorption curve). 

We use the prescription described above to compute the total 
energy absorbed by dust in the 'birth cloud' component (i.e. the 
dense central starburst) and in the surrounding 'ambient ISM'; 
this energy is then re-radiated by dust at infrared wavelengths 
(assuming conservation of energy). The total luminosity emitted 
by dust in the galaxy is then 



T tOt 



r BC j ISM 



(3) 



where L? c and L^ SM are the total luminosity re-radiated by dust 
in the birth cloud and in the ambient ISM components, respec- 
tively. 

This approach, d eveloped in ICharlot & Falll d2000t) and 
Ida Cunha et all d2008h . allows us to interpret the infrared emis- 
sion from galaxies consistently with the emission at shorter 
wavelengths using a simple energy balance argument (given that 
the main energy source in galaxies are stars, i.e. there is no con- 
tribution by an AGN). For normal star-forming galaxies with 
moderate amounts of dust (optically-thin case), the effective V- 
band optical depth ty seen by stars younger than tsc inside 
birth cloud is constrained by the attenuation suffered by emission 
lines, e.g. Ha/H/3, while the effective V-band optical depth seen 
by older stars in the diffuse ISM fit y is mainly constrained by 
the attenuated ultraviolet and optical continuum emission from 
the galaxy. 

In starburst-dominated ULIRGs, practically all the radia- 
tion from stars in the birth cloud component is absorbed by 
the optically-thick dust, and therefore it is not possible to reli- 
ably constrain ty using observations in the ultraviolet or opti- 
cal. In such cases, where the optical depths are very high, dust 
will absorb effectively even mid-infrared radiation, in particu- 
lar in the silicate band located around 9.7 fim, where the grain 
cross-section is maximum (and also to a lesser extent at 18 fim). 
The SpitzerflRS mid-infrared spectra of our ULIRGs allow us 
to measure the strength of the 9.7 fim silicate feature for each 
one of them. This feature strength can be converted to an appar- 
ent dust optical depth T9 .7, which can be used to compute the 
V-band optical depth by assuming dust optical properties and 
geometry. With this in mind, for each galaxy, we compute the 
strength of the 9.7 -fim silicate absorption feature, Ssi9.7//m, from 
the observed IRS spectra, as: 



> Si9.7//m 



In 



p obs 

17 cont 
r A 



(4) 



where F° bs is the observed flux density of the feature and F? ont 
that of the continuum, both evaluated at the wavelength, A$i, 
where the absorption presents its maximum (typically about 



Table 2. Properties derived from the IRS spectra of our ULIRGs. 
(a) strength of the silicate absorption feature; (b) V-band 
optical depth comp uted using the ULIRG silicate model of 
ISirockv etaf] f2008) assuming a cold dust screen; (c) equivalent 
width of the 6.2-ytmi PAH feature 1 . 
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9.7 fim). To compute the F^° nt we fit the IRS spectrum of 
each galaxy to a simple linear function with anchors at 5.5 and 
13.2 fim (rest-frame) and interpolate it through t he absorption 
feature . This method is similar to that used by Sirock vet alj 
(2008) for measuring the silicate strength in PAH-dominated 
IRS spectra, but here we use a long-wavelength anchor at 
13.2 fim instead of at 14.5 fim and a straight line instead of a 
power-law to fit the spectra. We note that both methods yield 
similar results for Ssi9.7/jm- We have checked that our values of 
the silicate absorption stren gth, shown i n Table [2] are consis- 
tent with those found by Hao et alj d2007l) for a larger sample of 
local ULIRGs and are significantly larger than the silicate ab- 
so rption strength of local starbursts of lower infrared luminosity 
of Br andl et"aD d2006l) . as expected from Fig. [2] 

We assume that the observed silicate feature of our ULIRGs 
originates in the birth cloud component (which can be approx- 
imated by an o ptically-thick dust sc reen), and use the ULIRG 
silicate model of Sirockv et al.l(l2008l) to estimate the V-band op- 
tical depth Ty 1 in the birth clouds from the apparent 9.1-fim op- 
tical depth inferred from the s ilicate absorption fe ature strength. 
The adopted model assumes a Mathis et al. ( 1977) grain size dis- 
tribution and a cold lOssenkopfetalJd 19921) dust screen, which 



i mply 



23.6 T9.7, where Tg- 



1 Si9.7jtan 



(see more details 



in Siroc kvet al.l2 008). We include a 30% uncertainty associated 
with Ty derived in this way to account for uncertainties related 
to the measurement of S si 9.7 fim (about 10%) and model assump- 
tions. Table [2] shows the silicate absorption feature strengths of 
our ULIR Gs, and the V-band optical depth predicted using the 
model of ISirockv et all d2008h . We will use these values to help 
constrain our model fy, as described in Section U.2.21 



1 We measure the 6.2-yum PAH feature equivalent width by, firstly, 
fitting the continuum emission close to the feature with a linear function 
anchored at 5.8 and 6.6 fim. Then, compute the feature emission by 
integrating the continuum-subtracted spectrum from 5.9 to 6.5 fim, and 
finally we divide this by the continuum emission interpolated at 6.2 fim. 
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3.1 .3. Dust emission 

We distribute L^ c and ij SM in wavele ngth between 3 and 
1000 fim using four main components (see Ida Cunha et all 2008 
for details): (i) emission from polycyclic aromatic hydrocarbons 
(PAHs; i.e. mid-infrared emission features), (ii) mid-infrared 
continuum emission from hot dust with temperatures in the 
range 130-250 K, (iii) emission from warm dust in thermal equi- 
librium with adjustable temperature, and (iv) emission from cold 
dust in thermal equilibrium with adjustable temperature. 

In stellar birth clouds, the relative contributions to L^ c by 
PAHs, the hot mid-infrared continuum and warm dust are kept 
as adjustable parameters. These clouds are assumed not to con- 
tain any cold dust. In the ambient ISM, the contribution to L^ SM 
by cold dust is kept as an adjustable parameter. The relative pro- 
portions of the other three components are fixed to the values 
reproducing the mid-infrared cirrus emission of t he Milky Way 
(more detail in section 2.2 of da Cunha et alj|2008l) . 

3.2. Median-likelihood estimates of physical parameters 

We use the model described in the previous section to inter- 
pret the observed spectral energy distributions of our ULIRGs 
and derive median-likelihood estimates of their physical param- 
eters. T o do so, we adopt a B ayesian approach similar to that 
used by |da Cunha et al.lp008l (see also jKauffmann et"al] 2003, 
iBrinchmann et al.ll2004lda Cunha eTallboiOl) . This method re- 
lies on a large library of random models encompassing all plausi- 
ble parameter combinations (e.g. star formation histories, metal- 
licities, dust optical depths, dust masses, dust temperatures) . 
Such a library was already built in Ida Cunha et al.l (12008, 2010) 
to derive statistical constraints on the physical parameters of lo- 
cal star-forming galaxies of moderate infrared luminosities. For 
each galaxy, we build the likelihood distribution of any physi- 
cal parameter by evaluating how well each model in the library 
accounts for the observed SED of the galaxy. The underlying as- 
sumption of our method is that the model library represents the 
distribution from which the data were randomly drawn, so it is 
important that the prior distributions of the parameters sample 
well the observational space and do not give too much weight to 
a priori implausible corners of the parameter space. 

To reproduce the observations of our ULIRGs, we find that 
some of the parameter priors have to be changed with respect 
to the priors used to model g alaxies with Lir < 10 12 L0 
in Ida Cunha etafl d2008l l2O10h . This confirms the idea that 
ULIRGs are probing a very different region of the parameter 
space of galaxies. In particular, we include higher dust opti- 
cal depths and temperatures and increase the probability and 
strengths of bursts of star formation. We describe the priors in 
detail in the next section. 

3.2.1. Model library 

We first build a large library of 50 000 stellar population mod- 
els with random star formation histories, metallicities and dust 
contents. We distribute the galaxy age uniformly between 2 and 
13.5 Gyr. We parametrize the star formation history of each 
model in terms of two components: an exponentially declining 
star formation rate of the form exp(-yf), where we distribute 
the time-scale p arameter 7 between and 1 in the same way as 
Ida Cunha et alJ (2008), and superimpose random bursts. Since 
our ULIRGs are highly star-forming galaxies, likely to be expe- 
riencing star formation bursts, we assign a 75 per cent probabil- 
ity of a burst occurring in the past 2 Gyr. The amplitude of each 



burst is defined as the ratio between the mass of stars formed 
during the burst and the total mass of stars formed by the un- 
derlying continuous model, A = Mb urs t/^cont- We distribute A 
logarithmically between 0.1 and 10. The duration of each burst, 
?burst> is distributed uniformly between 3 x 10 7 and 3 x 10 8 yr. 
The metallicity is uniformly distributed between 0.2 and 2 times 
solar. 

To account for the high optical depths of ULIRGs, we have 
to change significantly the priors of the attenuation parameters 
f v (effective V-band optical depth seen by stars younger than 
tsc in the birth cloud component) and fi (fraction of fy con- 
tributed by the diffuse ISM compone nt), compared to the pre- 
vious study of Ida Cunha et al. (2008). This is motivated by the 
very high infrared-to-optical ratios observed in the galaxies for 
which optical observations are available, and also by the high 
V-band optical depths t?, 1 derived from the IRS spectra using the 
silicate feature (Table [2}. The strength of the silicate feature in 
our ULIRGs implies that the birth cloud component has to be 
optically-thick, and the fact that ultraviolet and/or optical emis- 
sion is observed for several of our galaxies means that the dif- 
fuse ISM component should have much smaller optical depths. 
Therefore, we distribute fy uniformly in log between 2 and 50. 
These high optical depths imply that the optical emission from 
ULIRGs comes from a more diffuse, optically-thin component - 
the diffuse ISM component. An optically-thin ISM component is 
the only way the optical fluxes can be fitted consistently without 
the diffuse ISM contributing too much to the total infrared lumi- 
nosity (i.e. too high - Z/[ SM /L™). This is the motivation for 
drawing the V-band optical depth in the diffuse ISM, fify from a 
Gaussian prior centered at 1 with 0.5 dispersion, between and 
2. We note that this is in agree ment with studies of the ex tinction 
in spatially-resolved ULIRGs ( Garcia-Mann et al . 2009), where 
values of optical extinction as low as 0.2 mag are found for the 
external regions (in typical scales of 1 - 3 kpc). 

We also generate a library of 50 000 inf rared emission mod- 
els in a similar way to Ida Cunha e t al. (2008). We distribute the 
fraction of total dust luminosity contributed by the diffuse ISM, 
fn = Lj SM /L^ ot , uniformly between and 1. We randomly draw 
the fractional contributions by PAHs, hot mid-infrared contin- 
uum and warm dust in thermal equilibrium to the total infrared 
luminosity of b i rth clo uds, from the same prior distributions as 
Ida Cunha et al.l d2008l) . The equilibrium temperature of warm 
dust in birth clouds, T^ c , is distributed uniformly between 30 
and 60 K. For the equilibrium temperatur e of cold dust in t he dif- 
fuse ISM (7^ SM ), we extend the prior of Ida Cunha et all (2008) 
to include higher values in order to account for the observed hot- 
ter far-infrared colours of our ULIRGs; T^ SM is uniformly dis- 
tributed between 15 and 30 K. 

Finally, we combine the library of stellar population mod- 
els and the library of dust emission models by joining models 
with similar paramet er and scaling to the to tal infrared lumi- 
nosity L^ ot as detailed in lda Cunha et alJ d2008l) . For each model 
spectrum in our library, we compute the synthetic photometry in 
photometric bands from the ultraviolet to the far-infrared at the 
redshifts of our ULIRGs. 



3.2.2. Spectral fits 

We perform spectral fits by comparing the observed spectral en- 
ergy distributions and also the V-band optical depths of the birth 
cloud component derived from the silicate strengths measured 
using the Spitzer/IRS spectra in Section l3.1.2l to every model in 
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our stochastic library. Specifically, for each ULIRG, we compute 
the^- 2 goodness-of-fit of each model j in the library: 

where L' v and L' v , are the luminosities in the i th band of the ob- 
served galaxy and the / h model, respectively, cr, is the obser- 
vational uncertainty in L' v , w j is the model scaling factor that 
minimises^ 2 [see eq. (33) of Ida Cunha et al.l (12008)1. Ty is the 
V-band optical depth derived from the observed 9.7 fim silicate 
feature and cr si is the uncertainty associated with this measure. 
As discussed in Section [3. 1.21 we explicitly fit Ty to the effec- 
tive V-band optical depth in the birth clouds because they are so 
optically-thick this information is not accessible from ultraviolet 
and optical observations. 

We fit all the available observed broad-band fluxes of our 
galaxies from the ultraviolet to the far-infrared, plus the V-band 
optical depth of the birth clouds inferred from the mid-infrared 
silicate absorption strength (eq.|5]). We do not attempt at repro- 
ducing the detailed shape of the observed IRS spectra of our 
galaxies because our SED model is not designed to reproduce 
the detailed emission in this spectral range (so that the number 
of free parameters is kept minimal). However, we aim at repro- 
ducing the rough shape of the mid-infrared SED of our galaxies, 
since this reflects the emission by hot, stochastically-heated dust, 
and also by PAHs. To do so, we compute broad-band fluxes in 
the mid-infrared from the observed IRS spectra and include them 
as observables in the fitting procedure. 

We build the likelihood distribution of any given physical pa- 
rameter for our observed galaxies by weighting the value of that 
parameter by the probability exp(-^ 2 /2). Our best estimate of 
the parameter is the median of the resulting probability density 
function and our confidence interval the 16th-84th percentile 
range. The results of our SED fitting for the sample of 16 star- 
forming ULIRGs using this method are presented in the next 
section. 

4. UV-to-IR SEDs and the physical parameters of 
star-forming ULIRGs 

In Fig. [3] we show the results of our SED fitting for each of our 
16 ULIRGs. We compare, in the top panels, the best-fit mod- 
els (in black) to the observed multi-wavelength SEDs (in red). 
The residuals are shown in the bottom of each spectrum. The 
blue line shows the unattenuated emission from stars. For all 
our galaxies, this emission is very strong and blue, indicating 
that the stellar component is dominated by young stellar popula- 
tions resulting from current or very recent star formation. Most 
of this stellar radiation is absorbed by dust and re-emitted in the 
infrared, leading to the large observed infrared-to-optical ratios 
in these galaxies. The infrared emission of our galaxies peaks 
typically at wavelengths shorter than 100 /urn, indicating over- 
all high dust temperatures. For each galaxy, as a reference, we 
also plot the observed Spitzer/IRS spectrum (green line). We 
note that, even if we did not intend to fit the detailed shape of 
the mid-IR emission, our best-fit SED and the IRS spectrum are 
in good agreement in most cases. Both the shape of the PAH 
emission and the slope of the hot mid-infrared continuum are 
reproduced remarkably well considering the simplicity of the 
model. However, the fits do not reproduce the observed silicate 
absorption feature because dust self-absorption is not included 




Fig. 4. Distribution of the difference between the V-band opti- 
cal depth inferred from the observed silicate strength Ty and 
the best-fit model V-band optical depth of the birth cloud com- 
ponent fy est , in units of the measurement error <xsj. We assume 
a 30% uncertainty to account for measurement errors from the 
IRS spectra and also for uncertainties in converting from silicate 
strength to V-band optical depth. 



in our infrared SEDs. This should have a negligible contribution 
in terms of the energy balance performed here. 

Even if we do not explicitly include the silicate absorption in 
our infrared SEDs, the silicate optical depth measured from the 
observed IRS spectra provides valuable information about the 
optical depth, as discussed in Section [3] We find that in heavily 
dust-enshrouded systems as our ULIRGs, constraining well the 
dust optical depth is fundamental in our SED fitting and has a 
high impact in the statistical constraints of star formation rates, 
stellar masses, and contribution by the diffuse ISM and the birth 
clouds components to the dust heating. To illustrate how well 
the V-band optical depth in the birth clouds is recovered by our 
best-fit model, in Fig.|H we plot the distribution of the difference 
between the observational (t v ') and the best-fit (fy) V-band op- 
tical depth for the galaxies in our sample. In most cases, the V- 
band optical depth is recovered within 2cr (i.e. 60% of r^, 1 ). This 
is a satisfactory result considering the uncertainties in deriving 
Ty from the observed IRS spectra. 

An advantage of our method is that it enables us not only 
to find a best fit model, but also to compute the full likeli- 
hood distribution of each model parameter and identify possi- 
ble degeneracies between parameters. This allows us to deter- 
mine how well the parameters are constrained and how the avail- 
able set of observatio ns affects these constr aints (see more de- 
tails in section 3.2 of Ida Cunha et al]l2008l) . The bottom pan- 
els of Fig. [3] show the likelihood distribution of several physi- 
cal properties of our ULIRGs derived from the comparison of 
the observed SEDs with our stochastic library of models as de- 
scribed in Section l3.2.2l We list the median-likelihood estimates 
and confidence intervals of some important parameters for each 
galaxy in Table [3] The likelihood distributions of Fig. [3] show 
that all the parameters are reasonably well constrained for most 
of our galaxies. However, we can see some differences in how 
well some parameters are constrained from galaxy to galaxy, and 
this is primarily due to the different set of available photometry 
for each galaxy. Clearly, the more available data points, the bet- 
ter constrained are the parameters. 

It is important to note how essential it is to include ultraviolet 
to near-infrared observations, even if the emission from ULIRGs 
is largely dominated by longer wavelengths. UV to near-IR ob- 
servations probe stellar populations older than 10 8 yrs which 
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Table 3. Median-likelihood estimates and confidence ranges (computed as the 16th - 84th percentile range of the likelihood distri- 
bution) of some relevant parameters derived from the multi-wavelength fits to the observed SEDs of our ULIRGs: (a) fraction of 
total dust luminosity contributed by the diffuse ISM; (b) specific star formation rate; (c) stellar mass; (d) dust luminosity; (e) dust 
mass; (f) equilibrium temperature of warm dust in birth clouds; (g) contribution of warm dust in thermal equilibrium to the total 
dust luminosity; (h) equilibrium temperature of cold dust in the diffuse ISM; (i) contribution of cold dust in thermal equilibrium to 
the total dust luminosity. 
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heat the diffuse ISM (i.e. stellar population which do not belong 
to the central, heavily obscured starburst). These observations 
help constrain the relative amount of dust luminosity heated by 
these stars ( f^), the dust attenuation in the diffuse ISM (fit y), the 
stellar mass (M„) and, as a consequence, the specific star forma- 
tion rate (^/ s ). This is visible in the case of IRASF10156+3705, 
for which no observations short-ward of 3.6 /im are available, 
leading to typically more extended likelihood distributions (and 
hence wider confidence ranges) for these parameters. We note, 
however, that even in this case, the star formation rate is still rel- 
atively well-constrained because, given the input optical depth 
derived from the IRS spectrum, the model determines what is 
the necessary amount of star formation to power the very large 
infrared luminosity of the galaxy (constrained by the infrared 
observations). 

In the infrared, given the wealth of available data at least 
up to 160 fim, the total dust luminosity (Ll 01 ) and the contribu- 
tion to the total emission by PAHs (£p^' H ), the hot mid-IR con- 
tinuum (^mJr), and warm dust in thermal equilibrium as 
well as the temperature of this dust (T^ c ), are well constrained 
in most of the cases. Cold dust emits in the far-infrared to sub- 
millimetre, and hence is constrained by the available observa- 
tions from 160 to 850 fim. In general, the contribution of this 
component to the total emission (t;}? 1 ) is well constrained by 
these observations. However, the equilibrium temperature of the 
cold dust component in our model (T^ M ) is relatively hard to 
constrain for our ULIRGs even with far-infrared data, because 
the contribution of this dust to the total infrared emission is small 
and hard to isolate from the dominating warm dust component. 



We note that this contributes to the broadening of the likelihood 
distributions of the total dust mass, since cold dust, even while 
contributing little to the total infrared luminosity, contributes sig- 
nificantly to the dust mass. 

We have compared our median-likelihood estimates of the 
total dust luminosity ZJj ot with the infrared luminosities empiri- 
cally derived using the IRAS fluxes, Lir, presented in Table Q] 
The two luminosities agree within typically 0.07 dex (with 
no systematic offset) for most of the galaxies. For only four 
galaxies, the difference between the two luminosities is bigger: 
IRASF08208+321 1, IRASF10156+3705, IRAS 17208-0014 and 
IRAS 19458+0944, although even in these cases the maxi- 
mum offset between the two luminosities is at most 0.14 dex 
(i.e. 38%). 

It is important to check how our stellar mass estimates com- 
pare with the results of previous studie s, specially considering 
the ver y large optical depths of ULIRGs. Rodriguez Zaurin etalj 
d2010 ) derived the stellar masses of 36 local ULIRGs using op- 
tical spectroscopy and found that ULI RGs have mod erate stel- 
lar masses < 1.4 x 10 1 ' M (using a lSalpeterlfl955l IMF and 
Bruzu al & Chariot 2003 models spectral synthes is models; see 
also, e.g. iGenzel et al.ll2.Q01b iTacconi et alj|2002f). We n o te tha t 
the use of a Salpeter IMF and the iBruzual & Chariol (2003) 
models can make the M» estimates systematically higher by 
0.15-0.3 dex compared to the use of a Chabrier IMF and the 
latest version of the Bruzual & Chari ot models used in our study 
(lBruzualll2007tlda Cunha et al.ll2010h . Nevertheless, we find that 
our estimates of the stellar ma sses of our ULIRGs are consis- 
tent with the previous result of Rodri guez Zaurin et al .1 (1201 Oh . 
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Fig. 3. Top panels: Best-fit models of the ULIRGs in our sample (observed-frame). The IRAS name of the galaxy is indicated in 
the top left-hand corner. The blue line shows the unattenuated stellar spectrum, and the black line shows the total SED including 
the attenuated stellar spectrum and the emission from dust in the infrared. The red squares are the available observed broad-band 
luminosities, and the green line is the observed Spitzer/IRS spectrum of this galaxy. The residuals (L° bs - L™ od )/L° bs are shown at 
the bottom of the panel. Bottom panels: Likelihood distributions of physical quantities derived from our fit to the observed SED: 
fraction of the total IR luminosity contributed by dust in the diffuse ISM (f^); total effective V-band absorption optical depth of dust 
(f v); fraction of the total V-band absorption optical depth of the dust contributed by the diffuse ISM (//); effective V-band absorption 
optical depth of dust in the ambient ISM (/if y); specific star formation rate (if/s ); stellar mass (M») and dust mass (M^); total dust IR 
luminosity (Z^ ot ); global contributions by PAHs (^p^ H ), the hot mid-IR continuum (^\ R ) and warm dust to L temperature 
of warm dust in stellar birth clouds (T^ c ); contribution by cold dust to Lj°' (t;^), and temperature of cold dust in the diffuse ISM 



except for IRASF10156+3705, for which we find a significantly 
larger M», but this may be attributed to the more uncertain esti- 
mate of M» for this galaxy due to the lack of optical and near- 
infrared data (as mentioned above). 



5. Discussion 

5.1. Comparison between star formation rate and infrared 
luminosity 

To test the validity of our approach, we compare our star for- 
mation rate estimates with the ones obtained fro m the infrared 
luminosity using the formula of lKennicuttl dl998l) . This formula 
assumes that all the infrared emission from the galaxy is pow- 
ered by recent star formation, i.e. all the energy emitted by young 
stars in the galaxy is absorbed and re-radiated by dust, and no 
dust is heated by old stellar populations (optically-thick star- 
burst). In our model, this corresponds to having very small val- 



ues of the parameter, i.e. the contribution by the diffuse ISM 
(heated by stars older than 100 Myr) to the total dust luminosity 
V? 1 is very small (see Table[3]). Also, to power the huge infrared 
luminosities of ULIRGs with young stars alone, the effective at- 
tenuation in the birth clouds must be very high, as well as the star 
formation rates. Therefore we expect the star formation rates of 
our ULIRGs to be similar to the ones derived using the Kennicutt 
formula. 

In Fig.[5]and Table|U we compare our median-likelihood es- 
timates of the star formation rates of our ULIRGs with the star 
formati on rates pr e dicted from the total infrared luminosity us- 
ing the iKennicuti (0,998) formula for optically-thick starbursts 
(corrected for a Chabrier IMF). The star formation rate is com- 
puted as: 



M yr- 



1.08 x 10" 



(6) 
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Fig. 3. Continued. 
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Fig. 5. Median-likelihood estimate of the total dust luminosity 
plotted against star formation rate for the ULIRGs in our sam- 
ple (both derived from our SED fitting). The red line shows the 
Kennicutt (1998) relation for optically-thick starbursts, corrected 
for Chabrier (2003) IMF (eq. [6]). The error bars correspond to 
the 16th - 84th percentile range of the likelihood distributions of 
each parameter. 



where Lir is the total infrared luminosity between 8 and 
1000 fim, which we approximate with the total infrared lumi- 
nosity derived from our fits, L^ ot . 



We verify that our star formation rates agree remarkably well 
with the ones predicted using the infrared luminosities. The dif- 
ference between iff and ijj&z is typically less than 20% (except 
for IRASZ02376-0054 and IRAS17208-0014, but even in those 
cases i/<K98 is still in the confidence range of if/) . This confirms 
that the model is correctly attributing most of the infrared emis- 
sion of the ULIRGs to the starburst phase (i.e. the 'birth cloud' 
component in our model). This was only achieved through a cor- 
rect estimation of the optical depths enabled by the use of the 
IRS spectra. 

We caution that in more quiescent galaxies, dust heated by 
old stars in the ambient ISM contributes more significantly to the 
total infrared lumino s ity, an d the star formation rates obtained 
using Ida Cunha et all (2008) may diffe r sign ificantly than those 
derived using the formula of lKenn icutt ( 1998). In such cases the 
star formation rates would be overestimated using the Kennicutt 
law. 



5.2. The star formation histories of ULIRGs 

The model described in Section l3~T1 allows us to explore possible 
star formation histories of galaxies, and the effect they have on 
observable properties. One interesting result of the SED fitting of 
our ULIRGs is that, for all 16 of them, the star formation history 
corresponding to the best-fit model (Fig. [3) shows a strong burst 
of star formation occurring at the present age, i.e. all the best- 
fit models are clearly ones with a current starburst (see example 
of Fig. [6]). We find that, on average, this burst of star formation 
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Table 4. Comparison between the median-likelihood estimates 
of the star formation rate derived from our fits to the SEDs (a) 
and the star formati on rate compute d from the dust luminosity 
using the formula o f lKennicuttf l998 (b1. 



Galaxy 


^/M Q yr 1 
(a) 


(fes/Ms yr 1 
(b) 


IRASOO 199-7426 


mtZ 


210 


IRAS01494-1845 


186+Jf 


189 


IRASZ02376-0054 


230+™ 


304 


IRAS04114-5117 


218+«2 


180 


IRAS06009-7716 


98^ 


108 


IRASF08208+3211 


303!^ 


301 


IRASF10156+3705 


339 +72 


414 


IRAS 10565+2448 


»5_ 34 


100 


IRAS 12 112+0305 


230+0.0 


231 


IRAS13120-5453 


150!°;° 


176 


IRAS 16334+4630 


280!^ 


288 


IRAS 17208-0014 


qi + 102 

71 -56 


191 


IRAS 19297-0406 


941+10 


251 


IRAS 19458+0944 


160 +16 


168 


IRAS20414-1651 




170 


IRAS22491-1808 
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Fig. 6. Example of best-fit model star formation history of one 
of our ULIRGs, IRAS20414-1651, with a current age 3.56 Gyr. 



is responsible for forming 30% of the measured stellar mass of 
our ULIRGs over the last 100 Myr (this corresponds to about 
2 x 10 10 M© formed during the burst). 

The observed spectral energy distributions of our ULIRGs 
can only be well fitted with a simple, exponentially-declining 
star formation history with superimposed random bursts (as we 
use in our stochastic library of models; see Section 13.2. Il l if 
the star formation history contains a strong burst of star for- 
mation at the present age. This star formation history, while 
not necessarily being a unique solution, is a plausible scenario 
consistent with the main result that can be inferred from the 
SED fitting - that the galaxy must be experiencing a strong 
burst of star formation where a significant fraction of its stel- 
lar mass is being formed. To constrain the detailed star for- 



mation history of a galaxy, we would require optical spectro- 
scopic dat a to characterize the stellar popu lations of different 
ages (e.g. lRodrfguez Zaurin et al.l2008ll2010 ). Nevertheless, our 
results show that fitting the broad-band SEDs can at least pro- 
vide some clues on the star formation mode of these galaxies. In 
particular, our fits clearly show that, in order to reproduce the ob- 
served properties of our ULIRGs, a star formation history with a 
strong current burst of star formation is required. 

5.3. Comparison with local galaxies of infrared luminosities 
between 10 10 and 10 12 L 

One of the main advantages of the method used in this paper 
is that our results are directly comparable to results found for 
other samples of galaxies using the same framework. The phys- 
ical properties of a sample of 3258 local star-forming SDSS 
galaxies (z < 0.1) detected with IRAS, for which GALEX and 
2MASS observations are als o available (SDSS-/&4 S sample 
from now on), were derived in Ida Cunha et alj (1201 Oh using the 
technique used this paper. In this section, we compare the phys- 
ical properties of local galaxies in two bins of infrared luminos- 
ity: 10 < log(L IR /L Q ) < 12 (1538 galaxies with high S/N obser- 
vation s selected from the SDSS-IRAS sample of Ida Cunha et al.l 
2010), and 12 < log(L IR /L ) < 12.5 (15 ULIRGs from the 
present study). This is done in order to examine whether our 
findings based on SED modelling are consistent with the current 
understanding of ULIRGs. In Fig. [7] we plot the distributions of 
the median-likelihood estimates of several physical parameters 
related to the star formation activity and the dust properties of 
the galaxies. We describe the main similarities and differences 
between the two samples, even considering the low number of 
ULIRGs compared to the number of lower-luminosity galaxies. 

Fig- Gib) shows that the stellar masses of galaxies in the two 
infrared luminosity bins are similar, with the ULIRGs tending to 
be slightly more massive. We note, however, that the stellar mass 
estimates of ULIRGs are typically more uncertain than those of 
galaxies of lower infrared luminosities (and lower dust attenu- 
ation). In Fig. |7]c) we show that ULIRGs have typically much 
higher specific star formation rates than the SDSS-ZRAS galax- 
ies, implying doubling timescales of their stellar mass (assuming 
that they continue forming stars at the present rate) of the order 
of less than 1 Gyr [Fig. |7Jd)]. This is another confirmation of 
the well-known starbu rsting nature of our local ULIRGs (see, 
e.g. lGenzel etaDl998l) . 

We find that galaxies with 12 < log(L IR /L Q ) < 12.5 
also present very different properties of the dust emission than 
galaxies with 10 < log(LrR/Lg) < 12. Fig. [TJf) shows that, 
in ULIRGs, the contribution of the diffuse ISM component to 
the total infrared luminosity, f M , is typically very small (about 
0. 10), compared with w hat is found for the SDSS-IRAS galaxies. 
Ida Cunha etai1(l2010h find that, for these galaxies, dust heated 
by old stars in the diffuse ISM contributes significantly to the to- 
tal infrared luminosity, with most values of typically between 
0.40 and 0.70. 

Furthermore, Figs.|3g) and (i) show that the typical equilib- 
rium temperatures of dust grains in the 'birth cloud' component, 
r^ c , and in the 'diffuse ISM' component, T^ M , are typically 
higher in ULIRGs than in the lower infrared luminosity galaxies. 
Additionally, the contribution of warm dust to the total infrared 
luminosity fi? is very large for ULIRGs [Fig.Qh)], while the 
contribution of cold dust is much smaller than in the other 
sample [Fig. . Warm dust is the dominant contributor to the 
total dust emission in ULIRGs, and therefore, the overall tem- 
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Fig. 7. Comparison between the distributions of median-likelihood estimates of the physical parameters of 15 of our ULIRGs 
with 12 < log(ZJj ot /Lo) < 12.5 (red histograms), and of 1538 local star-forming galaxies of lower infrared luminosities (10 < 
log(L™/L Q ) < 12) selected from the SDSS-IRAS sample of Ida Cunha et al.ll2010l ( grey histograms), (a) dust luminosity, L™; (b) 
stellar mass, M»; (c) Specific star formation rate averaged over the last 100 Myr, if/$; (d) doubling timescale for the stellar mass, 
assuming the galaxies would continue forming stars at the present rate; (e) dust mass, Mj; (f) fraction of total dust luminosity 
contributed by the diffuse ISM, f^; (g) temperature of warm dust in the birth cloud component, T^ c ; (h) fraction of total dust 
luminosity contributed by warm dust in thermal equilibrium, fj?; (i) temperature of cold dust in the diffuse ISM, T^ SM ; (j) fraction 
of total dust luminosity contributed by cold dust in thermal equilibrium, 



perature of dust in our star-forming ULIRGs is higher than that 
of the SDSS-IRAS galaxies (see also lClements et alJl2010h . 

Based on the model fits we may also explore how the re- 
lations between the star formation rate, stellar mass and dust 
mass can be useful to study the star formation mode and evolu- 
tion of our galaxies. In Fig. [Ha), we plot our median-likelihood 



estimates of the star formation rate against the stellar mass for 
our 15 local ULIRGs with 12 < log(L IR /L ) < 12.5, and 
the 1538 SDSS-IRAS galaxies with 10 < log(L IR /L ) < 12. 
The SDSS-/RA5 galaxies present a strong correlation between 
ift and M„ (the Spearman test indicates a 18-cr significance level 
for this correlation). Previous studies of local and intermediate- 
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our ULIRGs do not show any statistically significant correlation 
between iff and M,. 

In Fig. [8fb), we plot the median-likelihood estimates of 
star formation rate a gainst dust mass for the two samples. 
Ida Cunha et alj d2010l) find a very tight correlation between if/ 
and M d for the SDSS-IRAS galaxies. Interestingly, our ULIRGs 
do not follow this relation. For a given dust mass, the ULIRGs 
have higher star formation rates than the SDSS galaxies - they 
are forming stars about 10 times more efficiently for a given 
Md. If we assume a con stant dust-to-gas ratio and rely o n the 
Schmidt-Kennicutt law dSchmidtl [19591: lKennicuttlll998h . this 
impli es denser star-forming gas in ULIRGs (as observed by 
e.g. iGao & SolomonlfcOOi iBrvant & Scovillel[l999t llono et al.1 
2009; I Juneau et al.ll2009h . These results are consistent with the 
scenario in which gas and dust in lo cal ULIRGs are c oncen- 
trated in a compact central region (e.g. Soifer et al. 2001; Diaz- 
Santos et al., 2010) following a merger of two gas-rich galaxies, 
and creating a strong starburst that powers the huge infrared lu- 
minosities of these systems ( as found by simulations of major 
mergers of spirals; e.g. iMihos & Hernquistlll994t lMihoslll999l; 
iHopkins et"atll2006l) . 

While these findings are not new, we show here how our SED 
modelling can be useful to derive relevant physical parameters 
which provide clues to the star formation mode of galaxies by 
using simple photometric observations. Additionally to a higher 
star formation efficiency, we also show dust in ULIRGs is typi- 
cally warmer, as shown in the histograms of Fig. [7] This is again 
consistent with the idea that the star formation in these galaxies 
is more spatially concentrated, and dust grains remain closer to 
the young, newly-forme d stars, reaching higher equilibrium tem- 
peratures (see also, e.g. lYang et alJ l2007; Gr oves et al.l 2008). 



Fig. 8. Median-likelihood estimates of star formation rate aver- 
aged over the last 10 8 yrs ft, plotted against stellar mass M, 
(a) and dust mass M d (b) for 15 of our ULIRGs with 12 < 
log(L™/LQ) < 12.5 (red points), and 1538 local star-forming 
galaxies of lower infrared luminosities (10 < l og(Lj ot /L(T)) < 
12) s elected from the SDSS-IRAS sample of Ida Cunha et alJ 
2010 (grey points). The errors bars correspond to the 16th-84th 
percentile range of the parameter likelihood distributions. The 
light grey lines show the fit to th e SDSS poi n ts: (a) best fit us- 
ing a slope of 0.77 as found in Elba z et al.1 (120071) : log(<A) = 
-7.467 + 0.77 log(M»); ( b) log (M d ) = 7.107 + LlllogOA), as 
found bv Ida Cunha et al.1 (120101) . 



redshift galaxies have found that star-forming galaxies exhibit a 
strong correlation between stellar mass and star formation rate, 
forming a tight 'main sequence' in which the range of pos- 
sible star formation rates of a galaxy is defined by its stellar 
mass (iBrinchmann et al.ll2004t iNoeske et al.ll2007bT: lElbaz et al.l 
120071) . A similar correlation has also been found for galaxies 
at z =! 1 (lElbaz et alj|2007b and at z 2 dDaddi et al.ll2007h . 
with the overall normalization of the correlation increasing with 
redshift, i.e. typical galaxies present high er star formation rates 
for a given stellar mass at high redshifts. INoeske et al.l ([2007a) 
explain this correlation in terms of a mass-dependent gradual 
decline of the star formation in disk galaxies as a result of gas 
exhaustion as they evolve. It is interesting to analyze the distribu- 
tion of local ULIRGs in the star formation vs. stellar mass plot. 
The ULIRGs are clearly above the relation defined by the SDSS- 
IRAS galaxies, i.e. at fixed stellar mass, ULIRGs have higher star 
formation rates (by about one order of magnitude). Additionally, 



6. Summary and Conclusions 

In this paper, we have studied the physical properties of a sample 
of 16 local star-forming ULIRGs observed in the mid-infrared 
with SpitzerflRS for which we have compiled additional pho- 
tometry from the ultraviolet to the infrared. 

We have interpreted the observed ultraviolet-to-infrared 
spectral energy distributions of these galaxies using the sim - 
ple and physically-motivated model of Ida Cunha et al.l (|2008). 
Using this prescription, we describe our ULIRGs using two main 
components: (i) an 'optically-thick starburst' with a lifetime of 
100 Myr, anal ogous to the 'bi r th clou d' component in the orig- 
inal model of Ida Cunha et a D J20081) . where dust is heated by 
young stars; (ii) a diffuse ISM component, with smaller dust ef- 
fective optical depths, where dust is heated by stars older than 
100 Myr. The model accounts for the emission of stellar pop- 
ulations from the ultraviolet to the near-infrared and for the at- 
tenuation by dust in the optically-thick starburst component and 
in the diffuse ISM. The infrared emission is simply computed 
by assuming that all the energy absorbed by dust in these com- 
ponents is re-radiated at mid- and far-infrared wavelengths, us- 
ing four main dust components. We have combined this model 
with a Bayesian technique to derive statistical constraints on pa- 
rameters such as the star formation rate, stellar mass, dust mass, 
dust temperature, and contribution of each dust component to the 
total infrared luminosity from their observed multi-wavelength 
spectral energy distributions. 

We show that the model of Ida Cunha et a D d2008t) is ver- 
satile enough to successfully account for the observed spectral 
energy distributions of star-forming ULIRGs. Although the ul- 
traviolet to near-infrared emission represents only a small frac- 
tion of the total power radiated by ULIRGs, observations in this 
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wavelength range are important to constrain the dust attenuation 
and the fraction of infrared luminosity contributed by the diffuse 
ISM, the total stellar mass, and the specific star formation rate. 
Additionally, our method can also provide clues to the star for- 
mation mode of our galaxies. The best fit models to the spectral 
energy distributions of our ULIRGs indicate that the observa- 
tions can only be optimally matched if the star formation history 
includes a strong burst of star formation at the present age. While 
this is not surprising, given the known starburst nature of these 
galaxies, it demonstrates the consistency of our method. 

To put the physical properties of local star-forming ULIRGs 
in context, we have compared them with a sample of SDSS 
galaxies of lower infrared luminosities (10 < log(LiR/Lro) < 12) 
detected in the infrared with IRAS (Ida Cunha et all 120 1 0) . Our 
results are consistent with a well-known fundamental difference 
in the star formation mode of local ULIRGs and other local 
galaxies of lower infrared luminosities. ULIRGs have star for- 
mation rates and dust temperatures significantly higher than the 
other galaxies, consistent with a merger-induced starburst, in 
which gas is driven to the central region, forming stars at very 
high rates and thus powering the huge infrared luminosities of 
these galaxies. 

The method used in this paper provides a framework to bet- 
ter understand local star-forming galaxies spanning a wide range 
of infrared luminosities in a statistical manner. This has the ad- 
ditional advantage that it can also be used to preform quantita- 
tive comparisons to high-redshift galaxies, for which often much 
more sparsely sampled SEDs are available. 

In future work, we plan to include the contribution by AGN 
in our modelling of the spectral energy distributions of galaxies, 
in order t o account for all types of LIRGs and ULIRGs in local 
samples dArmus et al.ll2009[ iDesai et al.ll2007l) . 

Acknowledgements. We thank the anonymous referee for comments and sugges- 
tions that significantly helped improve this paper. EdC, VC and TDS acknowl- 
edge partial support from the EU ToK grant 39965 and FP7-REGPOT 206469. 
We thank David Sanders and Vivian U for providing us the optical photome- 
try for some of our galaxies. We are grateful to Bernhard Brandl and Emanuele 
Daddi for useful discussions, and to Brent Groves, Nick Kylafis, Vivienne Wild 
and Stephanie Juneau for comments on the manuscript. This research has made 
use of the NASA/IPAC Extragalactic Database (NED) which is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 



References 

Armus, L., Charmandaris, V., Bernard-Salas, J., et al. 2007, ApJ, 656, 148 
Armus, L., Mazzarella, J. M., Evans, A. S., et al. 2009, PASP, 121, 559 
Blain, A. W., Smail, I., Ivison, R. J., Kneib, J.-P, & Frayer, D. T. 2002, Physics 
Reports, 369, 1 1 1 

Brandl, B. R., Bernard-Salas, J., Spoon, H. W. W., et al. 2006, ApJ, 653, 1129 

Brinchmann, J., Chariot, S., White, S. D. M., et al. 2004, MNRAS, 351, 1 151 

Bruzual, G. 2007, (astro-ph/0703052) 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Bryant, P. M. & Scoville, N. Z. 1996, ApJ, 457, 678 

Bryant, P. M. & Scoville, N. Z. 1999, AJ, 117, 2632 

Chabrier, G. 2003, PASP, 1 15, 763 

Chariot, S. & Fall, S. M. 2000, ApJ, 539, 718 

Charmandaris, V. 2008, in Astronomical Society of the Pacific Conference 
Series, Vol. 381, Infrared Diagnostics of Galaxy Evolution, ed. R.-R. Chary, 
H. I. Teplitz, & K. Sheth, 3-+ 
Chary, R. & Elbaz, D. 2001, ApJ, 556, 562 
Clements, D. L., Dunne, L., & Eales, S. 2010, MNRAS, 403, 274 
Clements, D. L., Sutherland, W. J., McMahon, R. G., & Saunders, W. 1996, 

MNRAS, 279, 477 
da Cunha, E., Chariot, S., & Elbaz, D. 2008, MNRAS, 388, 1595 
da Cunha, E., Eminian, C, Chariot, S., & Blaizot, J. 2010, MNRAS, 273 
Daddi, E., Dickinson, M., Morrison, G., et al. 2007, ApJ, 670, 156 
Desai, V, Armus, L., Spoon, H. W. W., et al. 2007, ApJ, 669, 810 
Dey, A., Soifer, B. T., Desai, V, et al. 2008, ApJ, 677, 943 



Downes, D. & Solomon, P. M. 1998, ApJ, 507, 615 
Due, P., Mirabel, I. F, & Maza, J. 1997, A&AS, 124, 533 
Dunne, L., Eales, S., Edmunds, M., et al. 2000, MNRAS, 315, 115 
Elbaz, D., Daddi, E., Le Borgne, D., et al. 2007, A&A, 468, 33 
Farrah, D., Afonso, J., Efstathiou, A., et al. 2003, MNRAS, 343, 585 
Farrah, D., Bernard-Salas, J., Spoon, H. W. W., et al. 2007, ApJ, 667, 149 
Gao, Y. & Solomon, P. M. 2004, ApJ, 606, 271 
Garci'a-Marfn, M., Colina, L., & Arribas, S. 2009, A&A, 505, 1017 
Genzel, R., Lutz, D., Sturm, E., et al. 1998, ApJ, 498, 579 
Genzel, R., Tacconi, L. J., Rigopoulou, D., Lutz, D., & Tecza, M. 2001, ApJ, 
563, 527 

Groves, B., Dopita, M. A., Sutherland, R. S., et al. 2008, ApJS, 176, 438 
Hao, L., Weedman, D. W., Spoon, H. W. W., et al. 2007, ApJ, 655, L77 
Hopkins, P. F, Hernquist, L., Cox, T. J., et al. 2006, ApJS, 163, 1 
Imanishi, M., Dudley, C. C, & Maloney, P. R. 2006, ApJ, 637, 1 14 
Iono, D., Wilson, C. D., Yun, M. S., et al. 2009, ApJ, 695, 1537 
Juneau, S., Narayanan, D. T., Moustakas, J., et al. 2009, ApJ, 707, 1217 
Kauffmann et al. 2003, MNRAS, 341, 54 
Kennicutt, Jr., R. C. 1998, ApJ, 498, 541 

Kennicutt, Jr., R. C, Armus, L., Bendo, G., et al. 2003, PASP, 115, 928 

Kim, D. & Sanders, D. B. 1998, ApJS, 119, 41 

Kim, D., Veilleux, S., & Sanders, D. B. 1998, ApJ, 508, 627 

Klaas, U., Haas, M., Miiller, S. A. H, et al. 2001, A&A, 379, 823 

Lahuis, F, Spoon, H. W. W., Tielens, A. G. G. M., et al. 2007, ApJ, 659, 296 

Laurent, O., Mirabel, I. F, Charmandaris, V, et al. 2000, A&A, 359, 887 

Le Floc'h, E., Papovich, C, Dole, H., & et al. 2005, ApJ, 632, 169 

Lutz, D., Spoon, H. W. W., Rigopoulou, D., Moorwood, A. F. M., & Genzel, R. 

1998, ApJ, 505, L103 
Lutz, D., Veilleux, S., & Genzel, R. 1999, ApJ, 517, L13 

Maddox, S. J., Efstathiou, G., Sutherland, W. J., & Loveday, J. 1990, MNRAS, 
243, 692 

Marshall, J. A., Herter, T. L., Armus, L., et al. 2007, ApJ, 670, 129 
Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 
Mihos, J. C. 1999, in IAU Symposium, Vol. 186, Galaxy Interactions at Low and 

High Redshift, ed. J. E. Barnes & D. B. Sanders, 205-+ 
Mihos, J. C. & Hernquist, L. 1994, ApJ, 431, L9 

Moshir, M. & et al. 1990, in IRAS Faint Source Catalogue, version 2.0 (1990), 
()-+ 

Murphy, Jr., T. W., Soifer, B. T, Matthews, K., & Armus, L. 2001, ApJ, 559, 
201 

Noeske, K. G., Faber, S. M., Weiner, B. J., et al. 2007a, ApJ, 660, L47 
Noeske, K. G., Weiner, B. J., Faber, S. M., et al. 2007b, ApJ, 660, L43 
Ossenkopf, V, Henning, T., & Mathis, J. S. 1992, A&A, 261, 567 
Rigopoulou, D., Lawrence, A., & Rowan-Robinson, M. 1996, MNRAS, 278, 
1049 

Rigopoulou et al. 1999, in ESA Special Publication, Vol. 427, The Universe as 

Seen by ISO, ed. P. Cox & M. Kessler, 833-+ 
Rodriguez Zaurin, J., Tadhunter, C. N, & Gonzalez Delgado, R. M. 2008, 

MNRAS, 384, 875 

Rodriguez Zaurin, J., Tadhunter, C. N, & Gonzalez Delgado, R. M. 2010, 

MNRAS, 403, 1317 
Sakamoto, K., Wang, J., Wiedner, M. C, et al. 2008, ApJ, 684, 957 
Salpeter, E. E. 1955, ApJ, 121, 161 

Sanders, D. B., Mazzarella, J. M., Kim, D., Surace, J. A., & Soifer, B. T. 2003, 

AJ, 126, 1607 
Sanders, D. B. & Mirabel, I. F. 1996, ARA&A, 34, 749 

Sanders, D. B., Soifer, B. T., Elias, J. H, Neugebauer, G., & Matthews, K. 1988, 

ApJ, 328, L35 
Schmidt, M. 1959, ApJ, 129, 243 

Scoville, N. Z., Evans, A. S., Thompson, R., et al. 2000, AJ, 119, 991 

Sirocky, M. M., Levenson, N. A., Elitzur, M., Spoon, H. W. W., & Armus, L. 

2008, ApJ, 678, 729 
Soifer, B. T, Neugebauer, G., Matthews, K., et al. 2001, AJ, 122, 1213 
Soto, K. T. & Martin, C. L. 2010, ApJ, 716, 332 

Spoon, H. W. W., Marshall, J. A., Houck, J. R., et al. 2007, ApJ, 654, L49 
Stanford, S. A., Stern, D., van Breugel, W., & De Breuck, C. 2000, ApJS, 131, 
185 

Shekel, M., Lemke, D., Klaas, U., Krause, O., & Egner, S. 2004, A&A, 422, 39 

Strauss, M. A., Huchra, J. P., Davis, M., et al. 1992, ApJS, 83, 29 

Surace, J. A. & Sanders, D. B. 2000, AJ, 120, 604 

Surace, J. A., Sanders, D. B., & Evans, A. S. 2000, ApJ, 529, 170 

Tacconi, L. J., Genzel, R., Lutz, D., et al. 2002, ApJ, 580, 73 

Vega, O., Clemens, M. S., Bressan, A., et al. 2008, A&A, 484, 631 

Veilleux, S., Kim, D., & Sanders, D. B. 2002, ApJS, 143, 315 

Veilleux, S., Rupke, D. S. N, Kim, D., et al. 2009, ApJS, 182, 628 

Yang, M., Greve, T. R., Dowell, C. D., & Borys, C. 2007, ApJ, 660, 1198 



Table A.l. Photometric data - ultraviolet and optical. All fluxes are in Jy. 
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Notes: (a) - SPSS Data Release 6 (model magnitudes); (b) - iMaddox etaP dl990T> : (c) - Howell et al. (2010); (d) -[S urace & Sanders (2000); (e) - D. Sanders and V. U, private communication; 
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Table A.2. Photometric data - near-infrared to mid-infrared. All fluxes are in Jy. 
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Table A.3. Photometric data - mid-infrared to far-infrared. All fluxes are in Jy. 
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4.68E-1 (A) 


8.18E0 W 


9.46E0 W 


9.30E-2 (,,) 


3.50E-1 (W 








2.81E-l ( s> 


4.90E-2 ( » 


IRAS 13 120-5453 


0.0307 


4.40E-1 W 


2.98E0 (C) 


4.11E1 (C) 


5.23E1 W 


6.42E-1® 


2.70E0' 6 ' 


3.13E1 W 










IRAS 16334+4630 


0.1909 


1.58E-2 (6) 


7.99E-2 (i) 


1.19E0 (n) 


2.09E0 (fl) 


2.10E-2 (,,) 


6.90E-2 (6) 


1.57E0 W 


9.40E-l (rf) 








IRAS 17208-0014 


0.043 


1.78E-1 (,,) 


1.61E0 (e) 


3.21El <e) 


3.61El <e) 


2.55E-1 < * ) 


1.21E0 (<,) 


1.91E1 W 






1.07E0* 


1.55E-1 W 


IRAS 19297-0406 


0.0857 


6.02E-2 (6) 


5.17E-l (i) 


7.32E0 W 


8.62E0 (<,) 


9.87E-2 (6) 


3.93E-1 W) 


6.64E0 (d) 


4.47E0 1 '" 








IRAS 19458+0944 


0.0999 


3.73E-2® 


2.53E-l (i) 


3.95E0 (,) 


7.11E0 (,) 


5A6E-2 (b) 


1.70E-l (£/) 


3.32E0 W 










IRAS20414-1651 


0.087 


2.51E-2 (6) 


2.64E-l (i) 


4.36E0 (n) 


5.25E0 (fl) 


3.89E-2 (6) 


2.20E-1' 6 ' 












IRAS22491-1808 


0.0772 


4.37E-2 (6) 


5.40E-l (f) 


5.54E0 W 


4.64E0 W 


7.90E-2 (6) 


4.09E-1 W) 


4.68E0 W) 


2.04E0 W 






1.90E-2" 



Notes: (a) - IRAS FSC dMoshir & et al. 1990); ( b) - Compu ted fro m IRS spectrum; (c ) - lYang et al.1 d2007l); (d) - Co mpute d using post - BCD products provided by the Spitzer Science Center, as 
described in lMarshall et al.ld2007h : (e) - [Sanders etafl d2003h : (/) - iDunne et all feOOCh : (g) - iRigopoulou et alJ d!996h : (h) - iKlaas et alj d2001h : (i) - IRAS P SC. 

Extra ISO photometry: IRASOO 199-7426 - F v (10//m) = 5.70F- 1, F v (120/zm) = 7.50 F0, F v (150/im) = 5 34£0, F v (180/zm) = 4.21E0, F v (200pm) = 3.00F0 dKlaas et al.l200lb ; IRASO 1494- 1845 
- F v (170jum) = 1.49F0 dStickel et al.ll2004h : IRAS 13 120-5453 - F F (170yum) = 3.05F1 dStickel et alj|2004h . 
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Appendix A: Multi-wavelength photometric data 



